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1. Introduction 

In this paper we analyze the pattern forming ability and pattern stability for a one-dimensional 
non-linear transport-diffusion equation on the circle. The distinguishing feature of this equation 
is the non-local turning velocity that is determined by interactions between particles in various 
orientations — velocity is given by a convolution term of an interaction rate V with the distribution 
function. In its general form, it also includes a diffusion term. 

In Section 2 we establish some basic facts on the transport-diffusion equation, like conservation of 
mass and symmetries, correspondence between solutions of higher periodicity for V and general 
solutions for its 'rolled- up' version V n . We also discuss the corresponding equation on the real line 
and its relation with the equation on the circle. Linearization near the constant stationary state 
provides conditions on the interaction rate V and on the smallness of the diffusion coefficient such 
that non-constant stationary states exist. These conditions are related to some conditions of Primi 
et al. [9], whose work deals with the same transport-diffusion equation. But Primi et al. succeeded 
in proving existence of n-peaks like steady states (for n > 1) using weaker conditions on V (so 
their result is stronger) using a method completely independent of linearization. Finally, we show 
that the constant stationary state is globally stable if the diffusion coefficient is large enough or 
the interaction rate small enough. For this proof we use the Fourier transformed version of the 
transport-diffusion equation. 

The method used by Primi et al. [9] to construct stationary solutions allows no conclusions on the 
stability of the pattern. Mogilner et al. [8] argued that without diffusion a single peak is stable 
if the interaction rate is everywhere attracting (this is a positivity condition on V). They used 
a discrete setting with peaks of equal masses. In Section 5.1 we generalize that method to peaks 
with different masses and analyze the stability of n > 1 equidistant peaks of equal mass by a 
linearization argument. 

The transport-diffusion equation is closely related to a non-linear integro-differential equation on 
S , in which the non-linearity comes from interactions between particles, and in which particles 
jump instantaneously from one orientation to another. We will comment on these relations in 
the discussion in Section 7, which concludes this paper. For the limit of exact alignment in this 
jump process the stability of a single peak has been analyzed by linearization near the peak, see 
Geigant [.:]. In Section 5.2 we use linearization of the transport equation (without diffusion) near 
n-peak solutions to see whether they are stable or not. It turns out that a single peak is stable if V 
is everywhere attracting and that two opposite peaks are stable up to changes in masses if V unites 
attracting and repelling features. Most interestingly, n > 3 equidistant peaks of equal masses are 
in general not stable. Under appropriate conditions on V the number of peaks is invariant, but 
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positions and masses of the peaks may change slightly. A technical difficulty is that solutions are 
invariant with respect to translations, therefore 'stability' always means stability up to translations. 

The integro-differential equation with exact alignment has also been studied by Kang et al. [6]. 
They find that solutions converge in the sense of distributions to peak solutions if the support of 
the starting function consists of disjoint, small intervals and if the interaction is attracting. In 
Section 4 we prove that the same holds true for the transport equation, but for the transport 
equation the assumptions on V are more local. 

An important fact in all proofs of convergence to peaks is invariance of the first moment or barycen- 
ter. However, the first moment is only invariant for the linearized equations or if the initial function 
has sufficiently small support. Example 5.2 in Section 5.1 shows that the first moment (if defined 
in a naive way) is in general not invariant. This is caused by the discontinuity of the function 
9 given by identifying S 1 minus a point with ] — i, h[- The same is true for the related integro- 
differential equation. 1 In Section 2.5 we show how barycenters can be defined locally and when 
they are invariant. 

Section 6 presents an algorithm for fast computation of solutions for the transport-diffusion equa- 
tion. It is based on the Fourier transform of the transport-diffusion equation. The same method 
has been used by Geigant and Stoll [3] for the integro-differential equation. A second algorithm, 
namely an iteration scheme, is used to find stationary solutions of the transport-diffusion equation. 
It is essentially the scheme for which Primi et al. [9] prove convergence for suitable turning rates 
V, and indeed, in our computations we observe that their assumptions on V are necessary for 
convergence. 

In Section 6.3 we present a number of examples obtained using our numerical algorithms. The first 
example shows the simultaneous bifurcation of first and second mode; near that bifurcation point a 
stable mixed- mode solution and a backward bifurcation exist. The second example shows a stable 
two-peaks like solution where the two peaks are not opposite. In the third example stable one-peak 
and two-peaks like solutions exist at the same parameter values. This contradicts a conjecture of 
Primi et al. [9] that a certain feature of V (basically its shape and the sign of the primitive) allows 
conclusions on the number of developing peaks, namely one versus two. In the last example we 
show that pattern formation may occur even if V is nowhere attracting. 

2. The non-linear transport equation with diffusion 

Let S 1 = M/Z be a circle of length 1. If we denote by p : M —> S 1 = K./Z the canonical projection, 
then we have associated maps p* from functions on S 1 to functions on R, where p*(f) = f op is 
the associated 1-periodic function on K, and from (sufficiently fast decaying) functions on K to 
functions on S , where 

zGR,pO)=0 

Definition 2.1. A closed interval I on S 1 is a closed connected subset that is not all of S 1 . Then 
I = p(I') for some closed interval I' = [a, b] C K (such that 6 — a < 1), and we write I = [p(a),p(b)] 
and call a — p(a) the lower end and j3 — p(b) the upper end of /. If h is a function on S , we write 

h(0)d6= / h(6)dd= / h(6)d8= / p*(h)(x)dx. 



Therefore we think that equation (26) in Kang et al. [6] and conclusions based on (26) are not correct. The 
authors have informed us that an erratum is in preparation. 
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If 6, ip £ I, we write 9 - ip £ R for the difference 9' - ip' where 9', ip' £ I' are such that p{9') = 9, 
p(iP') = ip. 

If y : S 1 —> M. is a function and / = ]o, b[ C R is an interval such that p(I) ^ S l , we will (for 
simplicity) say that ( V > on ]a, b[ if > on p(I) (equivalently p*{V) > on J); similarly for 
half-open or closed intervals. In the same way, we write V(a) for V(p(a)) if a £ R. 

2.1. The equation on the circle. 

We want to model a process that describes the change of orientation of filaments over time. The 
orientation is given by an 'angle' e 5 1 . The density of filaments at time t > with orientation 
9 £ S* 1 is given by f(t,9). The filaments turn continuously where the velocity of turning is 
determined by interactions with other filaments on the circle. At the same time there is random 
reorientation. This kind of dynamics is described by the following transport equation with diffusion. 

Bf 3 2 f 3 

(1) ~{t t 9) = D^(t,9) + -{{V* /(*,■)) • /(*, .))(<?), 

where D > is the diffusion coefficient and (V * f){9) = J sl V(9 — ip)f(tp) dtp is the convolution 
of V with / and gives the negative velocity of turning of filaments with orientation 9. 

We assume that the interaction function V : S 1 — > R is odd, because interactions with filaments 
on opposite sides of 9 must have similar consequences. In particular, V(0) = 0, i.e., there is no 
repulsion or attraction of filaments with the same orientation, and V{\) = 0, i.e., there is no 
interaction with filaments of opposite orientation. The sign of V(9) is important. If V(9) > for 
some interaction angle < 9 < 1 then the two filaments move towards each other, we call this 
'attracting'; if on the other hand V(9) < for some < 9 < I then the distance between the 
filaments becomes greater, they are 'repelling each other'. For odd V E C°°(S 1 ) Primi et al. [ 
prove a-priori estimates by which unique existence of smooth solutions of equation (1) can be 
shown. 

The following easy statement will be useful. 
Lemma 2.2. Let VJe C(5 1 ) with V odd. Then 

(V*f)(9)f(9)d9 = 0. 

Js 1 

Proof. We have 

(V*f)(9)f(9)d9= f [ V(9-i>)f(^)dtPf(9)d9= f f V{9 - VO/WOT) # dS. 
is 1 Js 1 Js 1 Js 1 Js 1 

If we swap ip and 9 in the last integral, it changes sign (since V is odd); therefore it must be 

zero. □ 



The following proposition states some basic facts on equation (1). 

Proposition 2.3. Equation (1) preserves mass, non-negativity, axial symmetry with respect to 
any axis and periodicity of initial functions. Moreover, the solution space is invariant under the 
group 0(2) of translations and reflections on S . 

Proof. Preservation of mass and positivity are shown by Primi et al. [9]. The remaining statements 
follow from the observation that the operator on the right hand side of equation (1) is 0(2)- 
equivariant (for the reflections in 0(2), this uses that V is odd). □ 
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Since the partial differential equation ( ) lives on S , the equation turns into a discrete system of 
ODEs when it is Fourier transformed. We denote by 

fk = / f{0)e- 27rtk9 d9 for k e Z 

the /c-th Fourier coefficient of / : S 1 — > R. Since / is real, /& = /_&; if / is even or odd, then 
A- G K or /ft g iK, respectively. For differentiable functions one has (/')fc = ^ikfk, the Fourier 
coefficients of a convolution are (V * f)k = Vkfk, and the Fourier transform of a product is the 
convolution of the Fourier series, {f ■ g)k = X^ez fi9k-i- 

Hence the Fourier transform of the transport-diffusion equation (f ) is 
f k {t) = -(2nk) 2 Df k + 27TikJ2v l f l f k -i 

l£Z 

(2) =c k f k +A7rk v lflfk-l forfceZ, 
where the eigenvalues c k G K of the system and the ^eR are defined as 

(3) c fc = -(2irk) 2 D + 4irkf v k and u fc = / * V(9) sm(2nk9) d9 = \V k . 

Jo 

Mass conservation is reflected by the equation /o = 0. Using /_& = fk, the equations with k < 
are redundant. We have c,t = c„j, = cj e K because = — V-k- By scaling D and V one may 
assume that the mass is 1, i.e., /o = J sl f(t, 9) d9 = 1 for all t > 0, which we will do from now on. 



2.2. Stationary solutions. 

In this subsection, we collect some results on stationary solutions of equation (1). We begin with 
some estimates on the shape of a stationary solution. 

Proposition 2.4 (Estimates for stationary solutions). Any stationary solution of equation (f) 
with D > satisfies the following ordinary differential equation on S . 

(4) D fe = - {y * f) - f - 

Let C = maxV/D and assume that f > is a stationary solution of (1) with mass 1. Then for 
9i,02 £ S 1 we have 

f(0 2 ) < f(0i)e c ^- 6 ^ and f(9 2 ) > f(9 1 )e- c ^- e ^ . 

In particular, 

max f 

< e 



max / ^ a C/2 



min / 

and therefore max / < e c ^ 2 m\iif < e c ^ 2 and min / > e _c / 2 max/ > e^ c ^ 2 . 

In any maximum # max of a stationary solution we have 

dV maxU' 

m an?/ minimum 9 m i n of a stationary solution we have 

d 2 f min V 

^Q2 (^ min ) — fj /(^min)j 
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If D = in equation (1), then / is a stationary solution if and only if (V * /) • / = 0. 
These results can be interpreted as follows. 

i) If the diffusion coefficient D is large compared to V, then any stationary solution is near the 
constant solution (or only the constant solution exists). 

ii) The inequalities for max / have a large right hand side when D becomes small. That leads us 
to expect that with decreasing D solutions may become large and maxima may be sharp peaks 
(the curvature is large). 

iii) If the minimum of / is small, then it is wide (the curvature is small). 

Proof. (Sec also Primi ct al.([ ]) for the derivation of the ODE.) Let + ^((V * /) /) = 0. 

Then D% + {V * /)/ = c, where c is some arbitrary constant of integration. But the integral over 
S 1 of the first and second terms is zero (see Lemma 2.2 for the second term), so c = 0. 

If / is nonnegative, then 

^{6)<^maxV J fWdif>f(9) = Cf{6); 

Therefore by integrating from d\ to 62, we get 

/(0 2 )< /(0i)e c|e2 - ei1 ; 

the other inequality follows by symmetry. We can then take 6\ = m i n and #2 = $max to be points 
where / attains its minimum and maximum, respectively. 

Now, let max be a maximum of / Because f'(9 max ) — we have 

jjg((V * /) ' = (V * /)(0max)/(0max) + (V * /) (0 max )/'(0 max ) < (max V)f(9 m!ac ). 

Therefore, 

= Df"(6 max ) + dg{(V * /) • /)(0 max ) < Df"(6 max ) + (maxl/')M nax ). 
If #min is a minimum of /, then the estimate in # m j n can be proved in the same way. □ 

Primi ct al. [9] use the ODE (4) to set up an iterative procedure for approximating stationary 
solutions. See Section 6.2 below. 

We state a simple consequence of equation (4) . 

Proposition 2.5. Assume D > in equation (I). IfV is -^-periodic and odd, then any stationary 
solution of (I) must also be ^-periodic. 

Proof. V{6 + i) = V(0) for all 6 € S 1 implies 

(V * f)(6 + i) = / V{6+1- VO/W <# = / V(6 - VO/W dip — {V * f){6) 
Js 1 ' Js 1 

for functions / on S 1 . If / is a stationary solution of (1), then it is a solution of the ODE (4). 
Since D > 0, it follows that /'// = -£(V * /) is ^-periodic. This implies that /(0 + I) = 7 /(0) 
with some constant 7, and since / > 0, we must have 7 > 0. Since obviously 7" = 1, we have 
7 = 1, and / is ^-periodic. □ 
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2.3. Solutions with higher periodicity. 

Let n > 1 and V : S 1 — > R be odd and continuous. We will be interested in ^-periodic solutions 
of equation (1). To understand these, the following functions V n and V n will be useful. 

n— 1 . n— 1 ~ 

3=0 4=0 

and are continuous and odd; is ^-periodic. In particular, 

— ) = for < k < 2n. 
\2nJ 

The following result shows that instead of considering ^-periodic solutions of equation (1), we can 
consider solutions of (1) without higher periodicity, when we modify the parameters D and V 
accordingly. 

Proposition 2.6. Let n > 1 and V be odd. Then there is a 1-to-l correspondence between —- 
•periodic solutions f of equation (1) and solutions f of 

(5) ^(t, 9) = n 2 D^(t, 9) + -((K * f{t, •))/(*, ■))(*), 

namely via f(t,n9) — fit, 9). 



Proof. Let /(t, ■) be a ^-periodic solution of (1). Define f(t,9) = fit, -). Then / has mass 1, and 
we have 

n n—1 /.i±i n n ~ 1 r~ ft 

4=0 •'n ^o- 70 

= /">„(- -v)/m)#=- /V n (^)/(t,^)#=-(f n */(t,-))w. 

Jo ^Tl ' n Jo \ n / \ n/ n 

Also, g(t,0) = I§£ (*,£). Therefore, 

- " 2 ^0(^ ) + ^ * /(*■ ■))/(*. •)) w 

= " 2 ^||(t, (9) + ~ ((K * />, ■))/(«, ■)) (*)• 

The converse can be shown in the same way. □ 



This shows in particular that diffusion acts more strongly on solutions of higher periodicity. In fact, 
we have max7 n < nmax V, so the quotient C = thokV/D in Proposition 2.4 will be multiplied by 
a number < — . 

— n 

In terms of the Fourier transformed system (2), we have fk = for n \ k, fk = f„k, and (V n )k = 
nV n k- So we only look at the equations with k a multiple of n and replace nk by k to obtain the 
system corresponding to equation (5). 
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2.4. The equation on the real line. 

A similar PDE can also be considered with R instead of S 1 as the spatial domain, 
(6) f t (t, x)=D^(t,x) + ^ {{W * g(t, •)) • g(t, •)) (x) . 

Here W : R — >• R is odd and g(i, •) is assumed to decay sufficiently fast, so that the convolution 
W * g(t, •) is defined. Note that the convolution is here given by an integral over all of R. There 
is the following relation between equations (1) and (6). 

Proposition 2.7. Let V : S 1 —> R be odd and define W — p*(V) (which is just V considered 
as a 1-periodic function on R). Let g : [0, T[ x R — > R be a solution of equation (6). Then 
(t,9) i y f(t,0) — P*{g(t, -))(0) is a solution of equation (1). 



Proof. We have 



dt y,/ t-i at 

x:p(x)—9 



E (c0(*.*) + ^(o^*fl(*.o)(*)-fl(*,*))) 

a;:p(a;)— 

D^(t,9) + -^(J s V(9-i>)f(t,ij)di>- J2 9(t,x) 



x:p(x)—9 



Here we use that 



p*(V)(x-y)h(y)dy = / ^(p(x) - W # ■ □ 

i is 1 

The advantage of equation (6) over (1) is that it is easily shown to not only preserve mass, but also 
the first moment (or, equivalently, the barycenter) of g(t, •), whereas the notion of 'first moment' 
usually does not even make sense on S 1 . 

Proposition 2.8. Let g be a solution of (6). 

(1) For any a £ R, (t, x) t— > g(t, x — a) is again a solution of (6). 

(2) (t, x) i — y g(t,—x) is again a solution of (6). 

/>oo 

(3) / g(t, x) dx is constant. 



(4) / xg(t,x) dx is constant. 

J —oo 

Proof. The first statement is clear (the operator on the right hand side is equivariant with respect 
to translations) . Since W is assumed to be odd, the right hand side is also equivariant with respect 
to x i — V —x, which implies the second statement. For the third statement, we compute 

d,t 



g(t,x)dx= / — (D-£{t,x) + {W*g{t,-)){x)g{t,x))dx = Q, 
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using the decay properties of g. For the last statement, we have 
d f°° __ u N ^ f°° _ d (^dg 
dt 



xg(t,x)dx = I x—(D-^(t,x) + (W*g(t,-))(x)g(t,x) S )dx 
(^ffC*, x) + (W* g(t, •))(*)<?(*, xfj dx 
W(x - y)g(t, y)g(t, x) dy dx 



OO poo 



OO J —oo 



= 0, 

since W is odd, compare Lemma 2.2. □ 

Later, we will consider the case without diffusion (so with D = 0) in particular. In this situation, 
compact support is preserved. 

Lemma 2.9. Assume that W is bounded and that D = in (6). Let g > be a solution. If g(0, •) 
has support contained in [a,b], then the support of g(t,-) is contained in [a — Ct,b + Ct] for all 
t>0, where C = \\W\\ oo \\g(0,-)\\ 1 . 

Proof. Let g + (t,x) = g(t,x + b + Ct). We show that suppg+(i, •) C ]— oo,0]. This implies that 
supp g(t, •) C ]— oo, b + Ct]. The argument for the lower bound is similar. 



We have 



^(t,x) = d t g(t,x + b + Ct) + Cd x g(t,x + b + Ct) 
at 

= d x ((W * g(t, ■) + C) ■ g(t, •)) (x + b + Ct) 

= d x ((W + llWJIoo) * g(t, ■) ■ g(t, -))(x + b + Ct) 

= d x ((W+\\W\\ 00 )*g + (t,-)-g + (t,-))(x) 



Now 



d 

dt Jo 



g+(t,x)dx= I d x ((W+\\W\\ 00 )*g + (t,-)-g + (t,-))(x)dx 
o 

i(W +\\W\\ 00 ) * g+(t,-) ■ g + (t 7 ■))(()) 

{Wi-yi + WWW^g+fov) dyg + (t,0) 



< 0, 

since g+ > and W + \\W\\oo > 0. On the other hand, 

POO 

g+(0,x)dx = and / g+ (t, x) dx > , 



so we must have / g+{t, x) dx — for all t. □ 
If W(x) is positive for positive x, we can say more. 

Proposition 2.10. Assume that W is continuously differ entiable with W'(0) > 0, that W(x) > 
for x > and that D — in (6). Let g > be a solution such that suppg(0, •) C [a,b]. Then 
g(t, ■) has support contained in [a, 6] for all t > 0, and it converges to a delta distribution m5 c with 
m = 1 1 <7 ( , -)||i and mc — j xg(0,x) dx, in the sense that 

/>oo 

lim / h(x)g(t,x) dx = mh(c) 

t— >oo I _ 
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g(t,x)dx= / d x ((W*g(t,-))-g(t,-))(x)dx = -(W*f(t,-))(b)g(t,b) 

b Jb 



for all twice continuously differentiable functions h : K — > R. 

Proof. Without loss of generality, m — 1 and c = 0. We first prove the statement on the support 
of g(t, •). In a similar way as above in the proof of Lemma 2.9, we see that 

d 

dt 

So if g(t, x) > for some t > and a; > 6, we must have g(r, b) > and (W * g(r, -))(b) < for 
some < r < t. Let £o be the infimum of r > such that g(r,b) > 0. Then g(to,x) = for 
x > b, and it follows that (W * g(to, -))(b) > 0. By continuity, we will have (W * g(r, ■))(&) > 
and g(r,b) > for all sufficiently small r > i , so that the derivative above cannot be positive. 
So g(t, x) > for some i > and x > b is not possible. This shows that suppg(£, •) C ]— oo, b] for 
all t > 0. The argument for the lower bound is similar. 

We now consider the second moment of g(t, •). Let M{t) = x 2 g(t, x) dx. There is /i > such 
that xW(x) > fix 2 for all |x| < b — a. Then 
d f°° 

-M(t)= / x 2 d x ((W * g(t, ■)) • g(t, -))(x) dx 



-2 / x{W*g(t,-))(x)g(t,x)dx 

J —OO 

pOO pOO 

-2 / / xW(x - y)g(t,y)g(t,x) dy dx 

OO 

(x - y)W(x - y)g(t, y)g(t, x) dy dx 



— OO J —OO 
OO pOO 



oo J —oo 
b r b 



<-fij / (x-y) g(t,y)g(t 7 x)dydx 
= -2fi [ x 2 g(t, x) dx = -2fiM (t) . 

J a 

This implies that M(t) < e~ 2 ^M(0); in particular, M(t) -> as t oo. 

Now let h € C 2 (R). We can write h(x) = h(Q) + h'(0)x + h 2 {x) where h 2 (0) = h' 2 (0) = 0. This 
implies that there is some C > such that |ft.2(^)| < Cx 2 for x <G [a, b]. We then have 

/oo 
(h(Q) + h'(0)x + h 2 (x))g(t, x) dx - h(0) 
-oo 

b r b 

h 2 (x)g(t,x)dx < / \h 2 (x)\g(t,x)dx < CM(t), 

J a 

and this tends to zero as t — > oo. □ 
2.5. Local masses and barycenters. 

For equation (1), the mass J s± /(f, 9) d6 is still an invariant, but there is no reasonable definition of 
a 'first moment'. (For this, one would need a function F : S 1 — > K that satisfies F(9 + a) = F(9) + a 
for all 9 £ S 1 and a € Ml. Such a function obviously does not exist.) However, we can define a 
localized version of a first moment. 

Definition 2.11. Let / : S 1 — > K be continuous and nonnegative, and let / C S 1 be a closed 
interval. Let 7' = [a, 6] C R be an interval such that p(F) — I. We define the local mass 

m{I,f)= ( f(9)d9= { P *(f)(x)dxeR 
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and, if m(I, f) > 0, the local barycenter 

M(IJ)=p(— i— / x P *(f)(x)dx) el. 

The local barycenter does not depend on the choice of I' — any other choice has the form I' + k 
with k € Z, and then we find that 

b+k c-b c-b 

xp*(f)(x)dx = / {x + k)p*(f)(x)dx = / xp*(f)(x)dx + km(IJ), 



' a-\-k J a J a 

so that the expression under p(-) changes by an integer, and the result is unchanged. 

Lemma 2.12. Let f > be a solution of equation (1), and let I C S 1 be a closed interval. Then 
the local mass m(I, f(t,-)) is time-invariant, provided there is no flow across the boundary of I: if 
a, /3 £ S 1 are the endpoints of I , then we require 

(V * /(*, •))(«)/(*, a) = (V * f(t, •))(/?)/(*, P) = 

for all t > 0. 

If in addition, there is no interaction with parts of f outside of I , meaning that 

V(0 - tp)f(t, 9)f(t, ip) = if Oel and ip I, 
then the local barycenter M(J, f(t, •)) is also time-invariant. 

Proof. We have 

jm(I,f(t,-)) = J de((V*f(t,-))f(t,-))(9)dO 

= (V * /(*, -))(P)f(t, 0) - (V * f(t, .))(a)/(t, a) = 
and, writing / for / and m{I , f) — m{I, f(t, •)) for simplicity, 

m(I,f)j t M(I,f(t,-)) 

xd x ((V*f(t,-))f(t,-))(x) dx 

b(V * f(t, -))(b)f(t, b) - a(V * f(t, .))(o)/(t, a) - [ (V * /(t, ^ 
V{6-i,)f{t,^)f{%e)d^d6 

1 

v(e-i>)f(t > i>)f(t,d)dipde 

i j i 

-o, 

because V is odd, compare Lemma 2.2. □ 
2.6. Invariance of support. 

We consider equation (1) without diffusion on the circle (but what we say here is also valid for the 
equation on the real line). 



US 
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Lemma 2.13. Let A = I\ U . . . U /„ C S 1 be a disjoint union of closed intervals in S ; write 
Ij = [a.j,(3j\. Assume that for every continuous function h : S 1 — > R+, we have the implication 

h\ A = => (V * h)(ctj) > and (V * h)(/3j) < for all 1 < j < n. 

Let f : [0, oo[ x S 1 — > M+ be a solution of equation (1) with D = such that /(0, •)!/! = 0. Then 
f(t,-)\ A = 0for allt>0. 



Equivalently, if supp /(0, ■) C S 1 \ A, then supp f{t, ■) d S 1 \ A for all t > 0. 

Proof. For the given solution /, let $ : [0,oo[ x S 1 — >• 5 1 be the flow associated to — (V * /), 

<P(0,9)=9 and ^(t,0) = -(7 * f(t, •))(*(*, 9)) . 
Then it is readily checked that for all a, j3 € S , the integral 

M(a,/3) = / /M)d0 

J$(t,Q) 

is independent of t. Write ay(t) = <&(i, «_,■), f3j(t) = Q(t,(3j) and 

n 

A(t)=\J[a j (t),(3 j (t)}, 
i=i 

then we have 

/ f(t,6)d0 = O 

JA(t) 

for all t > 0. Now assume that /(£, is not identically zero for some t > 0. Then we must have 
that A <f_ A(t). Let to be the infimum of all t > such that A <f_ A(t). Then for some 1 < j < n, 
we must have Oj(to) — ctj and ^-(to) > 0, or /3j(io) = /?j and ^ff(to) < 0. But in the first case 

^f(to) = ^(*o, aj) = -(V * /(to, ■))(*(*»>, = ~(V * /(to, •))(%) < , 
since /(to, OU = 0, an d similarly in the second case -34 (to) > 0, leading to a contradiction. □ 

3. Stability of the constant solution 

Theorem 3.1 (Local stability of the constant solution). The constant function f{0) = 1 is a 
stationary solution of equation (1). The eigenvalues of the linearization around f are 

Cfe = 4irk(—TrDk + v^) for k G Z, 

1 

where Vk — J 2 V(9) sin(27rfc#) d9. Hence, the constant stationary solution is locally stable if Ck < 
for all k > 0. 

More precisely, assume that ^2i >0 l\vi\ < 00 (this is for example the case when V is C 1 and 

I2(I2_ 



piecewise C 2 ); then p — y sup (>0 - ^ 6 ^ vf < 00. Define 

ii/ii 2 = E^ 



fe>i 



Any initial function /(0, •) such that 

II v., - max fc >i fc(f fc - far-D) _ min fc >i(-c fc ) 
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converges to the constant 1 in the sense that 

||/(V)-l||oo->0 and ||fl?/(t.-)||oo->0 for all n > 1 

as t — > oo. 
Remarks 3.2. 

i) The diffusion term has a stabilizing effect on the constant stationary solution. 

ii) Since by assumption, l 2 vi is bounded, c k will be negative for k 3> 0. Therefore, higher 
modes tend to be linearly stable. 

hi) Because periodicity is preserved, instability of the fc-th mode, i.e., c k > 0, implies that 
something interesting happens in the subspace of ^-periodic solutions. In all examples 
known to us there exist non-constant ^-periodic stationary solutions. However, we cannot 
exclude the possibility of time-periodic ^-periodic solutions or chaos. But see Corollary 4.2 
below, which shows that time-periodic solutions are impossible when there is no diffusion. 

iv) Since 



||/|| 2 = £i|/ fc | 2 <i £ |/jJ 2 = ^||/-i||2<^||/-i||L 

k>l keZ\{0} 

the condition on /(0, •) is satisfied when 



,1/(0,0 -i|U< 

Proof. Formulas (2) and (3) show that the given Ck are the eigenvalues of the linearization (which 
is obtained by disregarding quadratic terms on the right hand side of the differential equation). 

Recall that we assume fo — 1. We scale time by a factor 4tt and set S = irD in the Fourier 
transformed system (2) to get 

(7) fk = k((-dk + v k )f k + v ifih- 

*ez\{o,fc} 

Because = f k (/ is real) and = — v k £ R, we see that 
1 d^l 



2dt^r 1 ^k 

fe>l k>l 



k>l 

£>* - k5)\f k \ 2 +RcQr£t, i / i / fc jU 

k>l k>l l>k 

+ Re^£ £ vifif k f k _i - £ £ vJJ k f k+ i 

k>ll<Kk k>H>l 

5> fe -fa*)|/ fe | 2 +Re( J2 

^k-\-Tnfk-\-Tnfkfi 



k>l k,m>l 

Note that, setting k ■<— k + I, we have 

Re(£ £ vtfjkh-i) =Re(£ vififkfk+i) =Re(£ v lflfkfk+l I , 

fe>l l<i<fc fe,i>l fe,i>l 
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justifying the last equality above. In the remaining sum, we have set I <— k + m. We can estimate 
the latter as follows. 



Re^ ^2 v k+mfk+mfkf 
k,m> 1 

< | Vk+mfk+mfkfi 
k,m> 1 

EH 71 i \ fk+m fk fm 

^km(k + m)v k+m — —= —= 

'k,m>l W 

* (E \w) (E| E ^^> k+m ^= A 

k>l k>l m>l v ^ v 

fc>l fe>l m>l m>l 

= (E^iM 2 ) 2 (E(E^-^))^yi/ ; i 



fc+m I 



fc>l 1>1 k=Q 



fc>l Z>1 

t>u v 7 fc>l fe>l 

In terms of ||/|| 2 = £k>i j|/a;| 2 , this means 

(9) ^ll/II^Cmaxfe^-^+HI/lDll/ll 2 . 
Let /(0, •) satisfy the given condition and set 

c = -(maxk(v k - kS) +p||/(0,-)||) > 0. 

&;> 1 

It follows that 

||/(i)|| 2 < 11/(0, Oll 2 ^ 
for i > 0. Since l\f k \ 2 < ||/|| 2 , this implies that 

< >/k\\f(0, -)\\e- ct for fc > 1 and t > 0. 

We need a lemma. 

Lemma 3.3. Assume that \fk(t)\ < Ck a e~ ct for all k > 1 and a/Z i > 0, where C > and 
a < j are constants. Then for any to > 0, there is a constant C > (depending on to) such that 
\fk(t)\ < C'k a ~ 2 e- Ct for all k > 1 and aZZ i > t . 

Proof of Lemma 3.3. The quadratic part of the right hand side of the differential equation (7) 
for f k is 

-Rfe = v ififk-i + E v ihh-k - E y ififk+i 

0<Kk l>k l>0 
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We estimate Rk M . 



\Rk(t)\<C 2 e- 2ct ( \vi\l a \l-k\ a + J2\vi\l a \l + H a 



l>0 



< C 2 k a e- 2ct J2 C i l \ v i\ < C 2 k a e- 2ct < C 2 k 



a e~ ct . 



l>0 



Here we use that J2i>o^\ Vl \ < 00 ana - that l a \l±k\ a < C\k a l for some constant C\ only depending 
on a. Write c' k — k(v k — kS) = c k /(ATr) < 0. Then we have 



and therefore 



The integral is bounded by 



fk - c'kfk — Rk 
fk(t) = e c *7fc(0) + / e c * ( *- T >i4(T) dr. 



C,k a 



e -ct _ e c' k t 



<C,fc 



a-2-ct 



for some constant C 3 (note that — c' k k 2 and that \c + c' k \ = \c' k \ — c > p\\f(0, -)|| 2 > 0). This 



gives 



\fk(t)\<(k 2 - a e^+^\f k (0)\ + C 3 )k a - 



2 e~ ct . 



Since c' k + c < — const, k 2 , the first summand in brackets is bounded uniformly in k > for 
t > to > 0. This finishes the proof of the lemma. □ 



Repeated application of the lemma then shows that, given N > and to > 0, there is a constant 
Cn > such that 

I /*(*)! < C N k~ N e- ct for all k > 1 and all t > t Q . 

This implies 

||/(-,0-l||o O <2^|/ fe (<)| = O(e- ct ). 

k>l 

Similarly, for any n > 1, we obtain 

\\d$f(;t)\\oo < 2(2n) n Y, k n \f k (t)\ = 0( e - ct ). □ 

k>l 

Corollary 3.4 (Conditions for pattern formation). 

a) Instability of the first mode: If V is positive on ] 0, \ L i. e. all filaments attract each other, 
then the first mode is unstable for sufficiently small diffusion coefficient D: C\ > for 
D< 1. 

b) Instability of the second mode: Assume that there exists 9q £ jOjif such that V(9) > 
on ]O,0 Q [, V(6 ) = and V{9) < on ]0 O ,|[- Moreover let (*) V{6) > V{\ - 9) for 
9 G ]min(#o, \ — $0), \ [■ Then the second mode is unstable for sufficiently small diffusion 
coefficient D, i.e. c 2 > for D <C 1. 



Proof. Statement a) follows from v± > and ci = 47r(ui — ttD); statement b) follows from Propo- 
sition 2.6 and part a) because V2 > on ]0, |[. □ 
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The second condition (*) for instability of the second mode can be interpreted in the following way. 
If there are already two peaks forming then attraction towards the nearer peak must be stronger 
than towards the second peak. Using exactly the assumptions of b), it has been shown by Primi 
et al. [9] that equation (1) has a ^-periodic stationary solution with two equally large and very 
high maxima if the diffusion coefficient D is small enough. 



Theorem 3.5 (Global stability of the constant solution). Assume that 

,2 



k(k 2 - 1) 



P' = J2 '~ V k < 00 

fe>0 

(this is the case when V is twice continuously differentiable, for example). If 

Vk < nD k — — for all k > 1. 
k 

then every nonnegative initial function /(O, ■) £ C(S 1 ) of mass 1 converges to the constant func- 
tion 1 — there is some c > such that 

||W(V)-l)L = 0«(e- ct ) 

for all n> 0. 

For given V, the assumption is satisfied whenever 

p' + kv k 



D > max - , 

k>l k 2 TT 



i.e., if diffusion is strong enough. 



Proof. Inequality (8) implies 

(10) \ l ll/l12 " (E ^f^ v l]\fl\ 2 + f^Hvk - k6))\\f\\ 2 < y + maxfc( Ufe - kS)) ||/|| 2 . 
l>i 

(with ||/|| 2 = J2k>i il/fcl 2 as before and using |/;| < 1). Now the proof proceeds as for Theorem 3.1, 
but using (10) instead of (9). □ 



4. Convergence to peak solutions in case of small initial support and no diffusion 



If there is no random turning, i.e., D = in (1), then sums of delta peaks can be stationary 
solutions. To make this precise, we have to define the right hand side of equation (1) for suitable 
distributions on the circle. 

The kind of distribution we are mostly interested in are (positive) measures, but it turns out that 
it is advantageous to use differentiable measures instead. The main reason for this is that the map 
S 1 — > I? (5' 1 ), ip i— > S^,, is continuous, but not differentiable (since the derivative at zero would 
have to be —S' ). As a map to 2? 1 (5' 1 ), it becomes differentiable, though. 

Let k > 0. The space 2? fc (S' 1 ) is the dual space of the space C fe (S' 1 ) of k times continuously differ- 
entiable functions on the circle S . The elements of I? (5' 1 ) are called measures, and the elements 
of I? 1 (S' 1 ) are called differentiable measures on fir. By standard theory (see for example [5]), 
r D k (S 1 ) can be identified with the subspace of V{S 1 ) (which is the dual of C°°(S 1 )) consisting of 
distributions of order < k — f £ D k (S 1 ) if and only if 

k 

there is C > such that \{f,h)\ < ||fc w) ||oo for aR h £ C°° (S 1 ) . 

3=0 
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In particular, we can consider P°(S' 1 ) as a subspace of r D 1 (S 1 ). P(S' 1 ) can be identified with the 
space of 1-periodic distributions on M, compare [ ]. 

A distribution / 6 I?' c (5 1 ) is non-negative if (/, h) > for every non-negative function h € C 00 (5 1 ). 
We write (i* 1 ) for the set of non-negative distributions in V k (S 1 ). Note that for / £ D^S 1 ) 
and test functions hi < /12 we have (f,h\) < (/, /12). 

The support supp / of / G V k (S 1 ) is the smallest closed subset of S" 1 outside of which / = 0. 

Let '0 € S" 1 - The deZia distribution 5^ e P^S 11 ) is defined by (d^^h) = h(ip) where ft. is a test 
function. 

The moss of a distribution / e is defined as J sl f(9) dd — (/, 1). 

The convolution of a distribution / € T> k {S 1 ) with a function V € C^S 4 ) is defined as 

(V*/)(0) = </,V(.-^)> 
for 8 £ S 1 . It is known that Vt/eC'fS 1 ) (see [5]). For example, (V * 5^){9) = V(0 - V) for 

Convergence /„ —> f in I?(S' 1 ) (or V k {S 1 )) as n — > 00 means that (/«, /i) — > (/, ft) as n — > 00 for 
all test functions ft, e C°°(5 1 ). 

4.1. No time-periodic solutions. 

We now prove a lemma that is inspired by [ ] . It will let us deduce that there are no time-periodic 
solutions when there is no diffusion. We define 

$(0)= f F(0)#; 
Jo 

this makes sense as a function on 5 1 , since J gl V(V0 = 0- Note that $ is an even function. 
Lemma 4.1. Let f(t, 9) > be a solution of equation (1) with D = 0, and define 

*(*)=/ / ${6-i>)f{t,i>)f{t,e)d^d6. 

Js 1 Js 1 

Then ^ < 0, with equality if and only if f is a stationary solution. 

Proof. We have, by symmetry and integration by parts, 

j t *{t) = J J $(9 - ^)/(t, ^)^(t, 0) dfl + ^ J $(9 - i/>) f (t, i/>)f(t, 9) dip d9 

= 2 J J *(6-1>mt,i/>)9t(t,9)dil>de 

= 2j i ( y t $(0 - ^/(t, ^) <w)%({v * /(*, .))(«)/(<, ^))<w 

= -2 J — ($ * /(<, •)) (9) (v * /(*, •)) (e)f(t, 9) d9 

= -2 [ ((V*f(t 1 -))(8)) 2 f(t,9)d9 
Js 1 

< 0. 

From this computation, we also see that ^rr(io) = only if (V * f(to, -))f(to, ■) = 0, which implies 
§j(io; = for all 9 6 S 1 , so /(to, ■) must be a stationary solution. □ 
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In [8], this 'potential' ^> is introduced in the case when / is a sum of point masses. In terms of the 
positions of these point masses, equation (1) (with D = 0) then turns into a gradient system, which 
shows that it will evolve towards an equilibrium. If V is positive on 1 0, | [, then the potential has 
its global minimum when all the masses are concentrated at the same point. So Mogilner et al. 
conclude that the system of point masses will converge to this equilibrium (which corresponds to 
a single delta peak) if sufficiently perturbed from its initial state. We use the continuous analogue 
to show that at least no time-periodic solutions can arise. 

Corollary 4.2. Let f{t, 9) > be a time-periodic solution of (1) with D = 0. Then f is in fact a 
stationary solution. 

Proof. We have f(T, •) = /(0, ■) for some T > 0. This implies that *(T) = *(0). Since by 
Lemma 4.1, < 0, we must have ^rf(t) = for all < t < T. By Lemma 4.1 again, this implies 
that / is stationary. □ 

This argument does not carry over to the case with diffusion. Given the equalizing effect of 
diffusion, it seems rather unlikely that time-periodic solutions appear when D increases from zero. 
However, we cannot rule out the possibility of some kind of diffusion-driven instability that might 
lead to a periodic solution. 

4.2. Convergence to peaks. 

Let V e C 1 ^ 1 ) be odd, and let / € V^S 1 ). Let h € C 00 ^ 1 ) be a test function. Then the 
transport term of (1) is defined as 

(ii) (^((y*f)f),h) = -(f,(y*f)hf). 

Proposition 4.3. Let V £ C 1 (S 1 ) be odd and consider equation (1) with D = 0. 

i) A single peak f = € D^^S 1 ) with ijj € S 1 is a stationary solution of (1) with mass 1. 

ii) Let V(6o) = for fixed < 0q < \. Two peaks with arbitrary masses and distance 9o are 
a stationary solution, i.e. f — miS^ + m20~ x p+g G 2?i. (5* 1 ) is a stationary solution of (1) 
for any i/igS 1 and m\, TO2 > 0. If mi + m,2 = 1 then f has mass 1. 

iii) n > 3 peaks with equal masses and equal distances are a stationary solution: For 1 < j < n 
let ipj € S 1 with ip j+ i - ipj = i (where ip n+1 = ipt). Then f = i ^™ =1 G V + (S 1 ) is a 
stationary solution of (1) with mass 1. 

Note that n > 3 peaks with different masses are in general no stationary solution contrary to the 
"degenerate" case n — 2 (where V(h) = holds always). The situation changes e.g. for n = 4 if 
V(j) = 0. Then again stationary solutions consisting of four peaks with distance | and possibly 
different masses occur. 

Proof. Let ft, <E C°°(S' 1 ) be a test function. 

i) Using (11) and V(0) = 0, we get (64, {V * 5^)h') = V(0)h'(0) = 0. 

ii) Let ^2 = ip + . Using (11) and V(0) = V(0 ) = V(-9 ) = 0, we get 

(miS^i, + m 2 S^ 2 , (V * {m\8^ + m 2 (5^ 2 )) /l ') 

= mi{5^, (V * (mi^ + m 2 S^ 2 ))ti) + m 2 (5^, 2 , (V * (mi 5^, + m 2 S^ 2 ))h') 
= mi ( mi V(0) + m 2 V(-e ))h'(i/j) + m 2 { mi V(9 Q ) + m 2 V{0))ti \^ 2 ) = 0. 



iii) This follows immediately from statement i) and Proposition 2.6. 



□ 
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In the following theorem we show that solutions converge to single peaks if the support of the 
initial function is sufficiently small. 

Theorem 4.4. Let V G C 1 ^ 1 ) be odd with V'(0) > and V > on }0,9 V [, where < 9 V < \. 
Assume that f > is a solution of (1) with D = such that supp/(0, •) C I where I C S 1 is a 
closed interval with Vol(J) < 9 V and such that m(J, /(0, •)) = 1. 

Then supp /(£,•) C L for all t > 0, and f(t, •) converges to the delta distribution 5m, where 
M = M (I, /(0, •)) is the local barycenter of /(0, •) on I . 

Proof. We lift I to an interval I' = [a, t]cl and let £ = & — a < 9 V . Then /(0, ■) = p*(.go) for a 
function g : K — > M with suppgo C I'. 

We consider equation (6), where we take W — p* (V) on [—£, £] and extend it to all of K in such a 
way that it is odd and satisfies W{x) > for x > (which is possible since p*{V) > on ]0,£]). 
Let g be the solution of equation (6) with D = such that g(0, •) = go. By Proposition 2.10, 
suppg(i, •) C I' for all t > 0. So the function * -))g{t, •) appearing on the right hand 
side of equation (6) will always be equal to (p*(V) * ))g{t 7 •) (since — y) = p*(V)(x — y) 
when x,y £ I'). This means that 5 will also be the solution of equation (6), if we use p*(V r ) 
instead of W. By Proposition 2.7, we then have /(<, •) = p*(g(t, •)) for all t > 0. In particular, 
supp/(t, •) C p(suppg(<, ■)) C p(I') — I- By Proposition 2.10, we also know that g(t, •) converges 
to §m', where M' = J R xg(0, x) dx, so f(t, ■) — p*(g(t, •)) will converge to 5m, since M = p(M'). □ 

Example 4.5. Initial growth of an already sharp peak may be also seen if < D -C 1 even if no 
single peak solution is expected, see Figure 1. In Figure 1 the interaction function is ^-periodic, 
V{6) — sin(47r6'). Only the second eigenvalue becomes positive for decreasing diffusion, and Primi's 
et al. [9] conditions for existence of a single peak are not satisfied. Therefore, a single peak solution 
is not expected for (1) with D > 0, but as shown in the figure a sharp peak grows initially. Indeed, 
we did not see the development of a second peak although we had the program run up to times 
larger than 140. 

Note that by Proposition 2.5 any stationary solution must be ^-periodic in this example (and such 
stationary solutions exist and are expected to be stable). Therefore the behavior described here 
must be an artifact of the numerics. We think that the explanation is that the time scale for the 
transition from one peak to two peaks should be roughly of the order of e x ^ D , so the rate of change 
would be of order e~ x l D , which is numerically zero if D is as small as in the example. 

The next corollary follows directly from Theorem 4.4 and Proposition 2.6. The assumptions on V n 
imply also that the n-th eigenvalue is positive (Corollary 3.4). 

Corollary 4.6. Let n > 1, < d v < and let V G C^S 1 ) be odd and such that V^(0) > 
and V n > on ]0,0 V [. Let f > be a solution of equation (1) with D — such that /(0, •) is 
--periodic with mass 1, and assume that there exists an interval L C S 1 with Vol(J) < 9 V such 
that supp/(0,-)cU;r o 1 ( / +s)- 

Then the solution f converges to n peaks of equal masses at equal distances: 

/(*,■) 3 J 4 ast^oc, where M= A/(/,/(0, ■))• 
n 3=0 

The following corollary states that several peaks at random distances may form if V is zero in a 
neighborhood of i. There is, however, a minimal distance between them. 
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V=sin(4"pi'theta) D=0.005 



Figure 1. V{9) — sin(47r#), D = 0.005, and /(0, •) is the stationary solution 
of (1) for V(9) = sin(27r#), D = 0.005. The numerical algorithm is described in 
Section 6 (Fourier based method with 61 Fourier coefficients). 

Corollary 4.7. Let V £ C 1 ^ 1 ) be odd and < 9 V < \ such that V > on }0,6 V [ and V(9) = 
for 9 V < \6\ < i. Let n > 1 and /(0, ■) £ C+iS 1 ), supp/(0, •) C U"=i Ij where Vol(i» < 6 V and 
dist(Jj-, Jfe) > 0„ /or all I < j ^ k < n. Define m 3 (t) = m(Ij,f(t, •)) and Mf(i) = M(Ij,f(t, •))• 

T/ien supp(/(i, •)) C (J ■ Ij for all t > and t/ie masses mj as well as the barycenters Mj are 
constant in t. The solution f converges to a sum of delta peaks: 



mj S M . 



3 = 1 



as t 



Proof. As long as the support of /(t, •) stays contained in the union of the Lj, the evolution of / 
on each of the Lj proceeds independently, since the part of / contained in the other intervals does 
not contribute to the right hand side of equation (1). But then Theorem 4.4 shows that the part 



that starts in L stays in Lj and converges to a delta peak as stated. 



□ 



In the following theorem we are interested in convergence to two peaks, but Proposition 2.6 cannot 
be used since /(0, •) is not necessarily ^-periodic. Neither is V ^-periodic in general. 

We first prove a lemma that allows us to show convergence to the delta-distribution if mass is 
constant and second moments converge to zero. 

Lemma 4.8. Let JcS 1 be a closed interval. Let f n € P^S* 1 ) with supp/„ C I and (/ n ,l) = 
771 > for all ?i > 1. Let M G /, and let qM ■ S 1 — > K be a C°° function satisfying qu{9) = (0 — M) 2 
for all 9 £ L . If (/„, q^j) ~ > asn-> oo, then f n — > mSm as n — > oo. 

Note that the same conclusion is valid when we only assume that qui (9) > c(9 — M) 2 for all 9 £ I 
with some c > 0. 

Proof. Let i M : S 1 K be a C°° function such that £ M {9) = 9 - M for 9 £ L, and define 
a n = (fm^A-i) and b n = {fmlAi) Then by the Cauchy-Schwarz inequality (applied to the inner 
product (g, h) i— > (f n ,gh) for functions g,h : I M, concretely with g = 1 and h = £m) we have 
c? n < mb n . Since b n — > 0, we must have a n — > as well. Let h £ C°°(S 1 ) be a test function. We 
can write 

h{9) = h(M) + ti{M)l M {6) + r(9)q M (9) 
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for 9el, with \r(6)\ < C = \ max/ \h"\. We then have 

| (/» - m8 M , h)\ = \ (f n ,h(M) + h'{M)£ M + rq M ) - mh(M) \ 
= \h'(M)(f n (x),e M ) + (fn(x),rq M )\ 
< \h'(M)\a n + Cb n 

-40 asmoo. □ 



The final positions Mq and Mi of the two peaks in the theorem below are obtained from the special 
case when V(p(x)) = cx with c > in an interval around zero and V is ^-periodic. In this case one 
gets equations — ~ c Mj{t) (with Mj(t) defined as in the proof below), so that Mj(t) —> 0, 

justifying the choice of Mj. 

Theorem 4.9. Let V G C 1 (S' 1 ) be odd with V'(Q) > and V'(^) > 0, and assume that there exist 
< 6»i < 2 < \ such that V > on ]O,0 X [ and V < on]0 2 , |[. 

Let f > be a solution of equation (1) with D — such that /(0, •) G C + (S' 1 ) with J gl /(0, 0) d9 = 1 
and supp/(0, ■) C Jo U ii where Io CI S 1 is a closed interval such that Vol(/o) < min{0i, | — 02 } 
and Ii = Iq + i. T/ien supp /(t, •) C 7q U Jx / or a ^ i > 0. 

Let L be a closed interval in S 1 containing Iq U Ji. Define the local masses rnj(t) = rn(Ij, f(t, •)), 
and let M(t) = M(I, f(t, •)) be the local barycenter on I. Then mo(t), mi(t) and M(t) are constant 
in time; we write mo, m\ and M for their values. Define 

M = M + \m 1 and M x = M - \m =M -\. 

Then f(t, ■) converges to a sum of two opposite peaks: 

f(t, •) -» m a S A - Io + mxS^ as t^oo. 



Proof. We first show that supp/(i, •) C 7 U h for all t > 0. Let I = [a, (3), then I x = [a',/3'] = 
[a+ §]; let e = Vol(7 ) = f3 - a < min{0i, § - 6 2 }. Let h : S 1 R+ with supp/i C I U Ji. 

Because 7 > on ]0, 0i[U ] §, 1 - 2 [, we see that 

(V *&)(£) = / V(p-i/>)h(il>)dil>+ / V( )M^)#>0- 



o<»<ei i<.<i_e 2 

Similarly, 

(V*h)(a)= V( a-ip )h(ip) dip + V( a - ip )h(i>) dtp < 0, 

-Si<«<o e 2 <.<i 

because V is negative on both intervals. In the same way, we get (V*h)(a') > and (V*h)(j3') < 0. 
By Lemma 2.13 it follows that supp/(i, •) C Io U 1% for all £ > 0. 

By Lemma 2.12, the local masses m,j{t) are then constant, and the same is true for M(t) (since 
/(*, •) = on S 1 \ I for all t > 0). We now define local first and second moments by 

M j{t) = I {0- Mj)f(t,6)d6 and m 2j {t)= [ {9 - M d f f(t,6) dO for j G {0, 1}. 
hi hi 

Note that the expression 9—Mj makes sense on Ij (even on I — we lift to a suitable interval in R and 
compute the difference there). The definitions imply that Mo(t) + Mi(t) = M — mo-Mo — miMj = 
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for all t > 0. Let TOa(t) = rri2,o(t) + m2,i(i)- We will show that rri2(t) — > as t — » oo. The time 
derivative of m2 is (after integration by parts) 

^F ( * ) = _2 (/ / + / J i )v(d-i>)(6-M )f(t,i;)f(t,6)di,d6 

-2(jT jf +^ J i )v(6-i>)(6-Mi)f(t,i>)f(t,6)dipd9. 

To estimate this, we observe that there is 6 > such that 

V (</>)<£ > 60 2 and F(0 + |)0 > 60 2 for all e [-e,e]. 

This is because V > on ]0,e] and on ] |, | + e] and because V(0) = V(|) = 0, V'(0) > and 
V'{\) > 0. We now bound the various integrals from below. For the first, we find 

V(6 -i/>)(0- M ) f{t, iP)f(t, 9) dip dO 

V(e-ip)(6-M )f(t,i>)f(t,e)dil;d6+ f f V(i;-6)(TP-M )f(t,ilj)f(t,9)dipd6 

I a J Iq J la J la 



V(9-ip)(9-ip)f(t,i>)f(t,9) d^d9 
>b I f {9-^) 2 f{t^)f{t,9)d^d9 



'la Jlo 

= b[ [ ((9~M )-(ib-Mo)) 2 f(t,4,)f(t,9)dijd9 

Jlo Jlo 

= 2b(m 2 ,o(t)m - M (t) 2 ). 
In the same way, we find for the fourth integral that 

V(9-iP)(9-M 1 )f(t,ij)f(t,9) d^d9>b(m 2 ,i(t)m 1 -M 1 (t) 2 ) 



To Jh 



Hi Jh 

The remaining two integrals are estimated together, as follows. 

V(9-ij)(9-M )f(t,i>)f(t,9)dijd9+ f f V{9-^){9-M 1 )f{t^)f{t,9)d^d9 
/,,.// •// •//. 

V(9 - $){{B - M ) - (0 - Mi))f(t, </>)/(*, 9) dipd9 
>b I f ((9~M )-(4;-M 1 )) 2 f(t^j)f(t,9)d4>d9 

Jlo Jh 

= b(m 2 ,o(t)mi - 2M (t)M 1 (t)+m 2A (t)m ). 
Adding up, we find that (recalling that M (t) + Mx(t) = 0) 

^(t) < -26((m2,o(*) + m 2 ,i(t))(m + ™i) - M o(<) 2 - 2M (t)M 1 (t) - M x (*) 2 ) = -2bm 2 {t). 

This shows that rri2(t) < e~ 2b 'm2(0), and since m,2(t) > 0, this implies ni2(t) — > as t —> oo. 
So the local second moments 7712,0 (t) & n d J7i2,i(i) tend to zero as well. Using Lemma 4.8 on the 
intervals Iq,Ii separately, it follows that for t — > 00 the solution converges to two peaks, 

f(t, •) -> m 5M + mi 6^ , 

where the distance between the peaks is Mq — Mi = h- □ 



STABILITY OF PEAK SOLUTIONS 



22 



5. Linear stability of peaks 



5.1. Stability of position. 

We start this section by using the 'peak ansatz' of Mogilner, Edelstein-Keshet and Ermentrout [8]. 
The initial distribution is a sum of n > 2 peaks at positions 9j(0) € S 1 and with masses rrij > 
where X)j=i m j = 1 (different masses are a generalization of [8]). The solution keeps this form, 
f(t, 9) = Y^j=i m j^8j(t)(9), an( l the positions 6j(t) satisfy the following system of ordinary differ- 
ential equations. 

d9 n_1 

(12) -^■{t) = -Y i m k V{9 i {t)-e k {t)) for j = 0,...,n-l. 

fc=0 

To see this, we write 6g. = 5 {- — 9j) — S(. — 9j) and plug f(t, ■) = J2jrrijS(. — 9j(t)) into the 
transport equation d t f = dg((V * /)/)• For the left hand side we get 

(13) f(V)= X>i (-%)) S%-8j(t)), 

3=0 V J 

and for the right hand side 

(n— 1 n— 1 
£ E m jmk V(ej(t) - 9 k (t)) S(. - 9 3 (t)) 
j=0 fc=0 

(14) = E E mjmkV&it) - 9 k {t)) 5'(. - 9,(1)). 

j=Q k=0 

Comparing (13) and (14) we deduce (12). 

The case n — 2 is interesting. Since V(0) = and V is odd, the system is 

9 = -m V{9 - 9 X ) 

0'i = -m 1 V(9 1 - 9 ) = rmViOo - 9 X ) 
hence (recall that m + m\ = 1) 

(15) < L(e Q -9 1 ) = -V{9 Q -9 1 ). 

— 7} 

Because 9j(t) — > 9j implies 5g.( t ) — > 5§., we may conclude the following. 

Example 5.1. Let < 9 V < § and V € C 1 ^ 1 ) odd with V > on ]0, 9 V [ and V < on ]9 V , \ [. 
If dist(0 o (O), 9i(0)) < 9 V , then 9 (t) - 9i(t) for t ->■ cx), i.e., the solution of (12) converges 
to a single peak; if dist(#o(0), #i(0)) > 9 V , then dist(0o(i)> 9i(t)) — * §j hence the solution of (12) 
converges to two opposite peaks. dist(#o,0i) = 9 V is an unstable stationary solution. 

Now we are in a good position to show that one gets into trouble when defining a 'first moment' 

i 

in the 'obvious' naive way by fjj p*(f)(t,x)x dx. The point is that this is (in general) not time- 

2 

invariant. 

Example 5.2. Let < £ < §, V odd with V(0) = 9 on [0, 4e] and 

/(0,-) = |£e (o) + |£fii(o) where O (O) = + e) and X (O) = p(| - 3e). 

Then 

/(t, •) -> i (<*(• - (-1 - e )) + 5{. - (i - e))) £ - (I - e)) as f oo. 
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To see this, note that 

i) !»b.(o) = -|y(0 o (O) -0i(O)) = -|V(4e) < and ^(0) > 0; hence, dist(0 o (f), 6i{t)) is 
decreasing in t = 0; 

ii) !(l + o (t)-fli(t)) ( = 5 -V(\ + e {t)-9i{t)) = -(l + o (t)-0i(i)) as long as dist(0 o (t), 0i(t)) < 
4e. 

Since dist(0 o (O), 0i(O)) = 4e, i) and ii) imply that dist(0 o (t), 0i(*)) = 1 + 9 (t) - 9i(t) -> as 

i — > oo; 

hi) ^((9 (*) + 0i (*)) ( = ) 0, therefore, o (i) + 0i(f) = -2e (mod 1) for all t > 0. These facts imply 
that 0o (t) -> -5 - e and 0i(£) -> | - e. 

The 'first moment' of the initial distribution is p*(f)(0, x)x dx = —s. The 'first moment' of 

the limit is o~(a; — (| —e))xdx = ^ — e. Therefore, this 'first moment' is not invariant. (In fact, 
it jumps by | when one of the two peaks moves through the point p{\) € S 1 .) 

We will now analyze the local stability of two selected stationary solutions, namely peaks in one 
place, i.e., 9j — 9q for all < j < n, and peaks with equal masses at equal distances, i.e., mj = — 
and 0o 6 S 1 , 9j = 0j_i + A for 1 < j < n. Obviously, both are stationary solutions of equation (12). 
The matrix of the linearization is 
(16) 

( - "E m k V'(9 -9 k ) m 1 V'(9 ~9 1 ) ... rn n V> (9 - 9 n -i) ^ 

m o V'(9 1 -0 Q ) - E m k V'{9 1 -9 k ) ... m„_ 1 F'(0 1 - 0„_x) 

A = k=o,k=£i 

n-l 

, m o F'(0 n _! - O ) m 1 V'(9 n _ 1 - 9 1 ) ... - E m fc F'(0„_x - fe ) 

\ fc=0,fc^n-l / 

In both cases A has a clear structure such that the eigenvalues can be calculated explicitly (re- 
member E TO fc = 1 f° r the first case; if dj — 0j+i = — and rrij = — , then A is a symmetric and 
cyclic matrix, because V is even). The eigenvalues are 

r if 3 = 

Xj = < — V'(0) if 1 < j < n, assuming that 9 k — 9q for all k 

I £ Eti 1 V-'(*)(-l + cos(27rf )) if 1 < j < n, assuming that 9 k = 9 Q + 4 for all fc 
One eigenvalue is zero, because of the translational invariance of the system. Note that Xj = X n -j. 
Theorem 5.3. 

a) A single peak is stable up to translation in the space of peak solutions if and only ifV'(0) 
is positive. 

b) Let n > 2; n peaks with equal masses and equal distances are stable up to translation in 
the space of n-peak solutions with equal masses if Efe=i ^'(^X - ^ + cos (2tt^)) < (*) 
for all 1 < j < n — 1 . 

For n = 2,3 a necessary and sufficient condition for (*) ZioM is V'(— ) > /or a// 
1 < j < tt, — 1 ; the condition is sufficient for all n > 2. 

Proof. If V(0) > and V'(£) > 0, respectively, then all eigenvalues except Ao are negative; this 
implies stability up to translation. We now assume that (*) holds. 
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If n = 2, then > Ai = §V'(§)(-1 + cos(tt)) = -V'{\), so V'{\) > 0. 

If n = 3, then > X 1 = i(F'(i)(-l + cos(?f )) + V"(f )(-l + cos(f )) = f V"(|)(-l + cos(?f )), 



3V V3A ■"■ i 3 // 1 ' \3/* 1 \ 3 /' 3 V3A 1 "V 3 

; , , , 3 l) = V'(l), S oV'(l) = V'(l) 



because V'(§) = V"(i - 1) - V'(i), so V'(§) = V"(§) > 0. □ 



Example 5.4. Primi et al. [9] consider examples with V(ff) = sign(a) sm(2iv6 + a sin(27T0)) and 
find that four-peak like solutions are not stable if a = ±1.2. For D = this is explained now, 
because 

V'{\) = sign(a) cos(27r| + a sin(27r-j)) (2n + a2-K cos(27r|)) = -2n sign(a) sin(a) < 
for all < \a\ < n. Note that V'(|) > is still necessary for (*) to hold when n = 4. 

5.2. Stability with respect to small perturbations. 

We consider the linear stability of the stationary solution f(t, •) = S in the space of differentiable 
measures on S 1 , T> 1 (S 1 ). Recall that this is the dual space of C 1 (S' 1 ) and can be identified with 
the subspace of distributions in T>(S 1 ) of order at most 1. Note that Sg is close to S in V 1 {S 1 ) 
when 9 is small (since (Sg — 8,h) = Oh' (9) for some 9 between and 9, so that \\Sg — S\\ v i < \9\). 

We formulate a lemma that we will need later. 

Lemma 5.5. Let L be a (time-independent) differential operator on S , and let L be another 
differential operator on S 1 such that 

(Lf,h) = (f,Lh) 

for f G V(S 1 ) and h € C 00 ^ 1 ). Let fit,-) G V(S 1 ) be a solution of the PDE d t f = Lf. 
Let H(t,-) e C^iS 1 ) be the solution of the PDE d t H = LH such that H(0,-) = h. Then 
(f(t,-),H(—t,-)) is constant. Ln particular, 

(f(t,-),h) = (/(0, •),#(*,-)>. 

Proof. We have 

j t (f(t, -),H(-t, •)> = (d t f(t, -),H(-t, •)> + (/(*, ■), -(dtH)(-t, •)) 
= (Lf(t, -),H(-t, •)> + (f{t, ■), -LH(-t, •)) 
= (/(*,.), LH(-t,-)-LH(-t,-)) = 0. 

Applying this with h = H(t, •) instead of h to obtain H, we have H(—t, ■) = h and 

(f(t,-),h) = (f(t,-),H(-t,-)) = (/(0, •),#((),•)) - (f(0,-),H(t,-)). a 

Since the solution space of our equation is invariant with respect to translations, no stationary 
solution can be absolutely linearly stable. In order to deal with this technical problem, we will 
consider perturbations that do not change the barycenter. 

In the following, we will always consider equation (1) with D = 0. If we set f(t, •) = S+ f(t, •) and 
linearize, we obtain the linear PDE 

(17) f (t, 8) = ^(Vf(t, .) + (V* f{t, -))8)i9) . 
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Theorem 5.6 (Linear stability of a single peak). Let V £ C 1 (S' 1 ) be odd and such that V > 
on ]0, \ [ and V'(0) > 0. Assume that f(t, ■) £ I? 1 (S' 1 ) is a solution of equation (17) such that 
supp/(0, ■) £ I with a closed interval £ I £ S 1 with p(h) !■ We can lift I uniquely to an 
interval V £ R with £ I' and p(I') = I. We assume that (/(0, ■), 1) = (/(0, ■),£) = where £ 
is a function on S 1 that satisfies £(p(x)) — x for x £ V . Then f(t, •) converges to zero as t — > oo 
inV 1 {S 1 ). 

Proof. Let h £ C^S 1 ). We have to show that (f(t, •), h) — > as t — > oo. We have 

= (^(Wr) + (F*/(V))<^) 
= (Vf(t,-) + (V*f(t,-))5,-ti) 
= -(f(t,-),Vh')-(V*f(t,-))(0)h'(0) 
= (f(t,.),-V(h'-h'(0))}. 

(The last equality uses that V is odd.) We note that if h is constant, then (f(t, ■), h) is constant in 
time and that iih — ci on /, then the same is true. Since (/(0, ■), 1) = (/(0, •), £) = 0, (f(t, -),h) = 
for such h. We can therefore restrict to functions h satisfying h(0) = h'(0) = 0. Let H(t, •) denote 
the (unique) solution of the initial value problem 

f r -y% • , 

Then we see by Lemma 5.5 that (f(t,-),h) — (f(0, ■), H(t, ■)). (In particular, this shows that 
equation (17) has a unique solution in r D 1 (S 1 ) under the given assumptions.) Let $ : K x S 1 — > S 1 
denote the flow associated to V, i.e., 

Then H(t,$(t,9)) = h(8), as can be readily checked. Equivalently, H(t,9) = h($(-t,8)). Now 
we claim that H(t,-)\i converges to zero in C 1 (I). For this, note first that for 6 £ I, we have 
&(—t, 9) — > as t — > oo uniformly in 9 (this is because is the unique attracting and p(0) the 
unique repelling fixed point of the flow $). So (t, -)|/||oo — > \h(0)\ = 0. Next, we observe that 

For large t, $(—t,8) will be uniformly close to zero, so —V'($>(—t,8)) will be uniformly negative 
(recall that V'(0) > 0). This shows that J^$(-t,8) tends to zero as t — ¥ oo, uniformly for 9 £ I. 
This in turn implies that 



dH . 
~d9 



d 



(t,9) = h f m-t,8))—^(-t 



39 

also tends to zero uniformly on / as t — > oo. So 

(/(*, •), h) - (/(0, •), H(t, ■)> ->• as t -> oo, 

and this means that f(t, •) — > in I? 1 (S' 1 ). More precisely, it follows that supp/(i, •) C $(— i, /), 
so that the support is contracted to {0}, whereas mass and first moment are always zero. □ 

It is certainly natural to consider perturbations that do not change the total mass (thinking 
of redistributing the mass on the circle). What about perturbations that do not preserve the 
barycenter? Consider a small perturbation g in P 1 (5' 1 ) with mass zero and (g,£) — M with 
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\M\ <C 1. Then 5 + g = Sm + (S — Sm + g), and S — Sm + 9 is still a small perturbation, but now 
of Sm ■ Assuming that I — M, the theorem above then predicts convergence to the shifted 

peak S M - 

In a way, we can see this from the proof. If we do not assume that M — (/(0, •), f) — 0, then (using 
test functions h with h(0) = 0, but not assuming h'(0) = 0) we find that 

f(t,-)$-M8'. 
This is in accordance with Sm — S ~ —MS 1 . 

It is perhaps also interesting to compare Theorem 4.4 with Theorem 5.6. The former shows that an 
initial distribution that is contained in an interval covering less than half of the circle will converge 
to a delta peak under equation (1) without diffusion. The latter shows that this peak is stable with 
respect to small perturbations that avoid an arbitrarily small neighborhood of the point opposite 
to the location of the peak. 

Proposition 2.6 yields the following generalization to n equally distanced peaks with equal masses. 

Corollary 5.7 (Stability of n peaks with respect to ^-periodic perturbations). Let n > 1, let 
V £ C 1 ^ 1 ) be odd and such that V n > on]0, ^ [ and V£(0) > 0. Assume that f(t, •) € V 1 {S 1 ) is 
an ^-periodic solution of equation (17) such that supp/(0, •) C Uj=o «) with a closed interval 
£ J C S 1 with p(±^) ^ /. We can lift I uniquely to an interval J' C 1 with £ I' and 
p(I') = I- We assume that (/(0, •), 1) = (/(0, ■),£) = where £ is a function on S 1 that satisfies 
£(p(x)) — x for x G I' and £{6) — for 9 £ Uj=i + «)■ Then f(t, ■) converges to zero as t oc 
inV 1 {S 1 ). 

We now want to derive a result similar to Theorem 5.6, but for two opposite peaks of not necessarily 
equal mass. We take this stationary solution to be fo = m_<5_i/4 + m+Si/4,. If we take f = fa + f 
in equation (1) with D = and linearize, we obtain 

(is) I (< ' 9) = M ( {m - v{e + ^ + m+Vi ° ~ 2 ))/(t < 9) 

+ m-(V * f(t, -))(-i)*-i/4(fl) + m+(V * f(t, ■))(!) 

We will define 

V(6) = m^V(9 + i) + m+y(6» - ±). 

Theorem 5.8 (Linear stability of two peaks). Let V £ C 1 (S' 1 ) be odd and such that V'(0) > and 
V'(h) > 0. With the notations m_ 7 m + and V from above, suppose that V has exactly four zeros 
on S] namely —j, 9q, | and 9\ (in counter-clockwise order). Assume that fit,-) £ I? 1 (S' 1 ) is a 
solution of equation (18) such that supp/(0, •) C 2_ U L + with closed intervals ±4 £ L± C S 1 such 
that 7_ C p(]6»i - l,e {) and L+ C p(]6» ,6»i[). Let £ £ ^(S 1 ) be such that 1(9) = 9 - (±\) for 
9 £ L±. We assume that m(i+,/(0, •)) = m(/_,/(0, •)) = (/(0, ■),£) = 0. Then f(t, •) converges to 
zero as t — > 00 in I? 1 (S' 1 ). 

Note that V has to have at least four zeros, since V' is positive at the two zeros at ±|. 

Note also that because of the translational invariance, any two peaks with distance | are stable 
under the assumptions of Theorem 5.8. 



Proof. We proceed in a similar way as in the proof of Theorem 5.6. We find that 
d_ 

Ft 



(/>, -),h) = -(f(t, -),m-V(- + \){h' -ti{-\)) + m + V(. - \){h' - fc'(I))) 
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So we let H (t, 6) be the solution of 

f m) . !>(§ m - f <*.-») - D(f M) - § 

with iT(0, ■) = h; then (/(i, •), ft) = (/(0, ■), i? (t, •)) by Lemma 5.5, and we have to figure out the 
longterm behavior of H . We see that a function h that is constant separately on 7_ and on 7+ is 
a stationary solution on I_ U I + and that the same is true when h is a multiple of So we can 
assume that h{\) = ft(-j) = m+ft'(|) + m_ft / (-|) = 0. We write H' for Then we have that 

j t H'(t, \) = -m-V'{\){H'{t, \) - H\t,-\)) 

and 

jH'{t,-\) = -m + V'{l){H'{t, -\) - H% i)). 

This shows that m + H'(t, \) + m-H'{t,-\) = for all t and that H'(t, \) - H'(t, -A) -> as 
t -> oo (recall that V'(|) > 0). So iJ'(t, |) ->• and H'(t,-\) 0. By arguments similar to 
those in the proof of Theorem 5.6 (note that in the present situation, the flow associated to V 
moves the values of h away from | and — \ and toward $o and we then see that H(t, •) — s- 
in C 1 (/_ U /+) and therefore (f(t, ■), h) — > as t — > oo. For a general test function h, we then find 
that 

H(t, ■) -> ft(i)x+ + K-1)X- + (m+ft'(i) + m_/i'(-i))£ on /_ U 7 +! 
where x± is a function in C°°(S' 1 ) that takes the value 1 on I± and the value on I T . This 
translates into 

f(t, •) 4 m(I+, /(0, OMi/4 + m(I_,/(0, -))*-i/4 - (/(0, O^X^+^i^ + n*_^ 1/4 ) = 0. □ 

The need for the three assumptions m(I+, /(0, •)) = m(/_, /(0, •)) = (/(0, •), £) = arises because 
two opposite peaks of arbitrary masses and arbitrary orientation form a stationary solution. If we 
have a perturbation that violates these assumptions (but does not change the total mass), say 

m(J+,/(0,-)) = /i, m(/_,/(0, •)) = -/! and (f(0,-),£)=M, 
then we can proceed as in the one-peak case. We adjust masses and orientation to obtain 

(m + + n)S 1/4i+M + (m_ - m)£-i/4+M 

as a stationary solution such that the resulting perturbation of this solution satisfies the assump- 
tions. 

If there is some < 6 V < | such that V(6 V ) = and V'(8 V ) > (and V'(0) > 0, of course), then we 
expect two peaks at a distance of V also to be a stable stationary solution, up to a redistribution 
of mass between the two peaks and reorientation that preserves the distance. This is indeed the 
case. 

Corollary 5.9. Let V £ C 1 (S' 1 ) be odd and such that V'(0) > 0, and assume that there is < 
6 V < h such that V(9 V ) = and V'(9 V ) > 0. Let m± > with m + + m_ = 1. and consider 
the stationary solution fo — m+5g v /2 + « of equation (1) with D = 0. Let V(9) = 

m—V(9 + %■) + m+y(0 — f) and suppose that V has exactly four zeros on S] namely — 9q, 
%■ and 0i («n counter-clockwise order). Then fo is linearly stable with respect to perturbations 
satisfying the conditions in Theorem 5.8 (with ±| replaced by 

Proof. The proof is virtually identical to the proof of Theorem 5.8, after replacing ±i by ±%. □ 
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We saw that single peaks are stable up to reorientation if V'(0) > 0. Two opposite peaks are stable 
up to redistribution of mass and reorientation preserving the distance under the assumptions of 
Theorem 5.8, which include V'(Q) > and y'(^) > 0. Now we prove that n > 3 equal peaks at 
equal distances are stable if V (— ) > for < j < n — 1, up to ?. 

Theorem 5.10. Let n > 3. Let V € C 1 ^ 1 ) be odd and such that V'(£) > for all 0<j<n and 
V n > on ] 0, y~ [. Let, for < j < n. Lj be a closed interval in S 1 contained in ] ^ — ^ + ^ L 
and Zet £ G C°°(S' 1 ) &e sucft that £{9) = 9 - ^ for £ Ij, for all j. Then 

_^ n— 1 

tfte stationary solution fo = — <5j/„ o/ equation (1) uraift Z? = 

n i=o 

«s linearly stable with respect to perturbations f £ I? 1 (S' 1 ) sucft. iftai 

supp/ C M Ij, m(Ij,f) = for all < j <n, and (f,t) = 0. 

3=0 

Proof. The proof proceeds in a way analogous to the proofs of Theorems 5.6 and 5.8. The equation 
governing the development of H(t, •) is (writing again H' for ^) 

f)H 1 n 

(19) ^rfe f) = -~E ^ " ^( H '& ^ - £))■ 

The flow associated to = £^ ■ V(- — ^) moves away from the points ^ toward the points ^ ± 
So for any test function ft, satisfying = ^'(«) = for all j, we find that (/(i, •), ft) — > as 

£ — y oo in the same way as before. For the derivatives H'(t, -) we obtain the equation (using 

k=0 

This leads to 

j=0 j,fc=0 

with equality only if all £f'(£, ^) are equal. On the other hand, one sees easily that J2j H'{t, ^) 
is constant. Together, this implies that all H'(t, ^) converge to the same value as t — > oo. Since 
functions that are constant on each Ij and also i are stationary under equation (19), we get that 

n— 1 . n— 1 - n— 1 

3=0 j=0 j=0 

where \j € C°°(S' 1 ) is a function that takes the value 1 on Ij and the value on all Ik with k ^ j. 
In terms of /, this reads 

n— 1 1 n— 1 

/(£, ■) 4 £ m(/ J , /(0, - (/(0, ■), " E S 'j/n = 0- □ 

J'=0 3=0 

As before, if M = (f,£) ^ 0, then we expect a reorientation by M in the positive direction. It 
is less clear what happens when mass is redistributed between the domains of attraction of the 
various peaks. The proof above would suggest that we simply end up with equidistant peaks of 
different masses, but this will in general no longer be a stationary solution. If wc consider the 
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system of ODEs (12) for moving peaks, then we see by the theorem on implicit functions that 
there is a unique stationary solution up to translation near fo with peaks of prescribed slightly 
different masses if the matrix in (16) with 8j = 2- and rrij = — only has one vanishing eigenvalue 
— it is the matrix obtained as the Jacobian with respect to the Oj of the map 

n-l 

(m , . . . ,m n -i;9 , ■ ■ ■ ,9 n -i) 1 — > (- V'mfcV^ - 8 k )) , 

V k=o /o< j<n 

and the zero eigenvalue corresponds to an overall translation. This condition will be satisfied when 
the positions of n equidistant peaks of the same mass are stable up to translation, since then all 
the relevant eigenvalues are negative. Note that the condition on V in Theorem 5.10 is sufficient 
to ensure this is the case, compare Theorem 5.3. 

In the special case that we have V(^) = for all < j < n the stationary solution of the 
system (12) will consist of equidistant peaks even when the masses are not equal. In this case, one 
can formulate a variant of Theorem 5.10 in analogy to Theorem 5.8 that shows that n equidistant 
peaks with different masses are linearly stable with respect to perturbations respecting the mass 
distribution and the overall orientation. 

6. Numerical algorithms and simulations 

6.1. Solving the transport-diffusion equation via the Fourier transformed system. 

In Section 6.3 we will calculate numerically solutions of the transport-diffusion equation (1) for 
randomly chosen as well as pre-structured initial distributions. 

By using the Fourier transform we convert the partial differential equation into an infinite (but 
discrete) system of ordinary differential equations; since large Fourier coefficients of a smooth 
function are small we can then restrict to a finite system which can be solved very efficiently. 

In the Section 2.1 we found that the Fourier transform of the transport-diffusion equation (1) is 
given by (compare (2)) 

fk = -(2n) 2 k 2 Df k + 2mkJ2Vlflfk-l 

leZ 

(20) = c k f k + ink v ififk-i for k £ Z >0 , 

l£Z,l=£0,k 

where the eigenvalues c k of the system (see (3)) and v k € K are 

4 1 

c fe = -(2n) 2 k 2 D + Ankv k and v k = j V{8) sin(27r/c6>) dO. 

Jo 

Mass conservation is reflected by /o = 0; we put /o = 1. Note that the number of positive 
eigenvalues is usually small (see the remarks after Theorem 3.1). 

In order to avoid the necessity to use very small timesteps (k 2 is large for higher modes) we 
multiply (20) by exp(— c k t) and define g k (t) = fk{t) exp(—c k t). Then we get 

(21) g k (t) = (/ fc (t) - c k f k {t))e~ Ckt = 4^ e ~ Cfe * £ w,/i(t)/ fc _i(t), 

1^0, k 

which we solve by a second-order scheme. 

The number n of equations is adapted dynamically, in the following way. We start with n Fourier 
coefficients of /, assuming that higher modes are zero; we calculate the right hand side of (21) for 
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1 < k < 2n and accept for 1 < k < n the resulting (t + At) as the new value for fk ■ If for some 
k > n the slope of is larger than some (small) error bound, then the number n of equations is 
increased to fc + 1. The additionally needed Fourier coefficients with n < k < k + 1 arc initialized 
as zero. 

A further advantage of this scheme is that higher periodicity of an initial function is preserved. 

6.2. Solving the stationary equation via iteration. 

We also programmed the iteration scheme which Primi et al. [9] used to prove existence of peak-like 
solutions. 

We start with an arbitrary function /(°) on S 1 with given mass (usually 1). E.g., /(°) may be 
the solution of D d ^ ss -{9) = -V(9)f (0 \9), which is expected to lie near the one-peak solution, if it 
exists (see Primi et al. [9] ; V — V * S « V * f if / is one-peak like) . 

Then we iterate 

jf(n+l) p f 

D— (6>) = -(V^f {n) )(8)f {n+1) (9) andrequire / f {n+1 \9)d9 = / (0) (6>)d(9, 

d0 J S i J s i 

so that is a function on S* 1 with the same mass as f(°\ 

//this sequence converges, then the limit is obviously a solution of (4), i.e., it is a stationary 
solution of (1). Primi et al. [9] give criteria for convergence; e.g., the assumptions V(0) > and 
J V(ip) dip > for 9 € ]0, |] imply the existence of one-peak like solutions if D is small. 

Our iteration program reliably finds stationary solutions with ^-periodicity if no stationary so- 
lutions with lower periodicity are present. Otherwise, it is better to use V n and n 2 D instead of 
V and D, compare Proposition 2.G. (Unfortunately, in our program numerical instabilities accu- 
mulate, so that ^-periodicity of for n > 1 is not preserved numerically, in contrast to the 
theoretical prediction.) 

6.3. Examples. 

In the following examples the stationary solutions were calculated with both algorithms (exceptions 
will be mentioned); their stability was checked with the Fourier based system. 

The first example is interesting because it shows a backward bifurcation and mixed mode solutions. 

Example 6.1. Let D > 0, 7 > and define 

V{6) = sin(27r(9) + 7 sin(47r(9). 

We use formula (3) for the eigenvalues and get 

1 

ri 1 
ci = -47r 2 £> + 4tt / (sin 2 (27r6l) + 7sin(47r6») sin(27r6>)) dO = tt (-4ttL> + 1) > D < — 

Jo 4tt 

and 

1 

= -I6n 2 D + 8tt / ' (sin(27T0) sin(4n9) + 7sin 2 (47r6»)) dO = 2t:(-8ttD + 7) > <^ D < —. 

Jo 87r 

For (£>, 7) € M + X R + all four combinations of c\ <> 0, c 2 <> are possible, see Figure 2. All 

other eigenvalues are negative. 

For all parameter values we have V'(0) > and J® V(tp) dip > for < 9 < |; therefore for 
very small diffusion coefficient one-peak like solutions exist (Primi et al. [9]) and at least for 
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c\ < 0, c 2 > 

j/ C 2 = 



y = 2 ci < 
c 2 < 



ci = 



V has single zero 
l 

2 ' ' V 2 



V > 



Figure 2. Regions of in/stability of first and second eigenvalue 



D = they are stable by Theorem 5.6. We find that V 2 {0) = 7sin(47r6>), therefore K,'(0) > and 
Jq V2(ip)dip > on ]0, |[, thus two-peaks like solutions exist for small enough diffusion (Primi 
et al. [9]). For 7 < J we have V'(0) > and V'(|) < 0, therefore two peaks are not stable, 
see Theorem 5.3; two-peaks like solutions are only stable in the space of ^-periodic solutions (if 
D is sufficiently small), see Proposition 2.6 and Theorem 5.6 for D = 0; note that c 2 > 0. For 
7 > |, V has a single simple zero and V'(h) > 0, so for D — two-peaks solutions are stable by 
Theorem 5.8. 

For D = ^ rj 0.0796 and 7 = 2, first and second eigenvalue of (1) are zero simultaneously. 
Therefore stationary solutions with two maxima of different height can be expected to exist, so- 
called 'mixed mode solutions' (Golubitsky and Schaeffer [4]); Figure 3 (top left figure) shows such 
solutions. 

Figure 3 shows typical stationary solutions and their stability for 7 = 2 and 7 = 4. The ^-periodic 
stationary solutions are unstable (7 = 2), or they become unstable with decreasing D (7 = 4). 
This suggests that in general, the stability result for two peaks at D = cannot be carried over 
to (very small) D > 0. However, solutions need much longer times at smaller D to move beyond 
states with two peaks of different height. E.g., for D of size of the order of 0.05, 7 = 2, and starting 
with a small perturbation of / = 1, we get two peaks of different heights in the first two time units, 
while convergence to the mixed-mode solution needs about 30 time units; for D ss 0.01 and 7 = 2 
as well as 7 = 4, these time scales change to one unit for initial pattern formation and several 
hundred units for convergence to one peak. 

In Figure 3 the stationary solutions were generated with the iteration method, and their stability 
was tested with the Fourier algorithm. The unstable ^-periodic steady states in Figure 3 are stable 
in the space of ^-periodic functions; they are also found with the Fourier-based algorithm if the 
simulation is started with a ^-periodic function. 

The second example shows non-trivial solutions when V'(0) and V'(^) are negative, and hence for 
zero diffusion coefficient neither one peak nor two peaks at distance | are stable by Theorem 5.3. 

Example 6.2. Let 

V{6) = sm(27T0) - sin(47r60. 
Only the first eigenvalue is positive for small enough diffusion coefficient. Figure 4 shows how 
the stationary solution is approached. As one expects according to Corollary 5.9 (which holds for 
D = 0), for small diffusion coefficient it consists of two peaks with distance 6 V (where V(9 V ) = 0). 
We checked numerically that indeed V'(0) = \ (V'(0) +V'(9 V ))> 0. 
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1 /2-periodic 

6.80 D=0.07 (c1= 0.4, c2= 1.5) stable 




1 /2-periodic 

7.4 D=0.1 (c1=-0.8, c2=9.3] stable 




Figure 3. Stationary solutions for V{6) = sin(27r6>) + 7sin(47r#) with 7 = 2 
(top) and 7 = 4 (bottom) (V'(0) > 0,V'{\) > 0,V(0.29) = and F(0.27) = 0, 
respectively), (top left) Stable mixed mode solutions; the distance between the 
maxima is | . Note the backward bifurcation — there is a stable stationary solution 
even though both eigenvalues C\ and C2 are negative! For D > 0.093 and 7 = 2 
we found no non-constant stationary solution, (top right) A ^-periodic stationary 
solution; it is unstable, but stable in the space of ^-periodic solutions. For D > 
there are no non-trivial ^-periodic stationary solutions, (bottom left) Mixed mode 
solution with distance | between the two maxima; for D = 0.15 we found no 
stationary mixed mode solution, (bottom right) ^-periodic solutions. In mixed 
mode solutions the larger maximum grows with decreasing D while the second 
maximum vanishes. In ^-periodic solutions the maxima grow with decreasing D. 
In all cases the peaks become narrower. 

This solution could be calculated only with the Fourier transformed system, since the iteration 
method does not converge — it runs into a two-cycle. 

The next example shows that one-peak and two-peaks like solutions are possible in the same 
model at the same parameter values. It is interesting that there exists a one-peak like solution 
although the first eigenvalue (of the linearization near the homogeneous solution) is negative for 
all parameter values. 



Example 6.3. Let 



V{9) = sin(47T0) + 7sin(67r6»). 
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Figure 4. Development of the stationary solution for V(9) — sm(2n9) — sm(4n0) 
started with f(0,Q) = 1 + 0.8 cos(27r6>). 9 v = 0.17 is the only non-trivial zero of 
V. The numerical computation started with 20 Fourier coefficients and ended 
with 140. 



For — | < 7 < | the turning rate V has a single zero on ]0, |[, V'(0) > and V'{\) > 0; also, 
V'(|) > for 7 > |. If 7 > 0, then for all 5 > there is an e > such that / Q 8 V(tjj) dtp > e for 
S < 9 < |. Theorem 7.3. in Primi et al. [9] shows that a one-peak like solution exists for small 
enough D. The second eigenvalue C2 is positive for D < the third eigenvalue C3 is positive for 

D < 12¥> Which is ~ - 013 fOT 7 = - 5 - 

Figure 5 shows stationary solutions (left side; calculated with the iteration scheme) and how they 
are approached in time (right side; Fourier based program). We see that one-peak and two-peaks 
like solutions are locally stable for small enough D and 7 = 0.5. The one-peak like solution develops 
when the initial distribution is sufficiently centered (compare Theorem 4.4), but in the simulations 
/(0, •) did not have compact support. The three-peaks solution is stable in the subspace of ~- 
periodic functions; it is unstable for D = 0.01 (data not shown); For D = 0.001 a solution with 
three peaks of different height, of which only two had distance ^, developed when the simulation 
was started with a perturbed three-peaks like function. Note that V(3) actually satisfies the 3-peaks 
stability conditions of Theorem 5.10. 

With small diffusion coefficient the 'typical' outcome of a simulation that is started with small 
deviations from / = 1 are one large and one small peak that are opposite. We suppose that for 
D > these become equally high for large times; the smaller D is, the more time will be needed 
for that. 



Caption for Figure 5. Top row (left): These stationary solutions for various D- values were 
calculated with the iteration algorithm; for D 0.02 the solutions look ^-periodic; (right): A 
one-peak like solution with a small second maximum develops fast when the simulation is started 
with a centered distribution; here we started with the stationary solution for V(9) — sin(27r#), 
D = 0.03. 

Second row (left): These ^-periodic stationary solutions for various D-values were calculated with 
the iteration algorithm with forced ^-periodicity, (right) Approximation in time of a ^-periodic 
stationary solution. The simulation was started with f(0,6) = 0sin(27r(9), cj> = 0.1 (start not 
shown) . Until t m 8 the distribution converges towards the homogeneous distribution (this is what 
it has to do — any initial distribution that is orthogonal to the modes occurring in V dies out), 
then triggered by some numerical noise, instability of the constant solution takes over, and the 




FIGURE 5. V{9) = sin(47T0) + 0.5sin(67r#). 
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Figure 6. V{6) = sin(47r6>) - 3 sin(27T0) (left), V{6) = sin(67r(9) - 4 sin(2vr0) 
(right) D = 0.01. Initial condition were one (left) and two (right) sharp peaks. So- 
lutions were calculated with the Fourier-based algorithm. The iteration algorithm 
does not converge for these V's (it runs into two-cycles); however, it converges 
to the shown stationary solutions if V^),2 and ^3^3 and |- and ^-periodicity are 
used, respectively. 

second mode begins to grow. 

Third row (left): ^-periodic solutions were calculated with the iteration scheme with forced |- 
periodicity. (right) The simulation was started with a small perturbation (Re/i(0) = 0.01, 
Re /2(f)) = 0.005) of the ^-periodic solution for D = 0.01. A distribution with three slightly 
different peaks develops very fast, where the distances between first/second and first /third peak 
are not |; then no further changes are discernible. 

Bottom row: 'Typical' result of a simulation, here started with the ^-periodic solution for D = 0.01 
which was perturbed by Re/i(0) = —0.01, Re/2(0) = —0.005. A distribution with two different 
maxima at distance | develops fast, then no further changes are discernible. 
End of caption for Figure 5 

The last example shows that a variety of behaviors is possible if V < on ] 0, h [. 

Example 6.4. We compare 

V (2 )(0) =sin(47T0) -7sin(27T0) (7 > 2) and V (3) (0) = sin(6vr0) - 7 sin(27r0) (7 > 3). 

In both cases V < on ]0, | [, V(0) < and V'(h) > 0. For V( 2 ) om Y the second eigenvalue is 
positive for D < for F(3) only the third eigenvalue is positive for D < ; all other eigenvalues 
are negative. We have V(2),2 = 2 sin(47T0) > on ] 0, | [ and VWa = 3 sin(67r#) > on ] 0, | [; all 
other V(j\ n are zero. 

Therefore we expect (at small enough diffusion coefficient D) for V( 2 ) stationary solutions with 
two equal maxima at distance | , and for V( 3 ) three equal maxima with distance | . These develop 
indeed, but the time scales are interesting, see Figure 6. Two, resp. three different maxima develop 
very quickly but at unexpected distances; development toward equal distances and heights can be a 
very slow process. The explanation is that for both V there are orbits of other stationary solutions 
when D = — for V12) two peaks whose masses add to 1; for Vr^) three peaks whose positions and 
masses satisfy (12), namely = J2k=i m kV(9j(t) — 8k(t)) for j = 1, 2,3 and mi + m 2 + m 3 = 1 
(9j e S 1 , rrtj > 0). 
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7. Discussion 

For the transport-diffusion equation (1) a wide variety of different patterns has been observed. 
Indeed, only a limited number could be shown in the last section. It emerges that it is nearly 
impossible to predict pattern formation only by knowing the shape of V; however, if one compiles 
information like the sign of the eigenvalues Ck, the zeros of V, the signs of V (-) and the shapes of 
the V n and of f Q V n (ip) dtp, then the picture becomes clearer. 

If there is no diffusion, then we know quite something about the stability or otherwise of peak 
solutions. Stability of n peaks shows itself often also in the 'short-time' behavior of solutions of the 
diffusion-transport equation at small diffusion. Therefore it is hard to clarify numerically whether 
a given stationary solution is stable for small diffusion. It is an open and interesting problem how 
to clarify the stability of stationary solutions if diffusion is present and if several eigenvalues are 
positive. 

A possible interest in the transport-diffusion equation (TDE) comes from its relation to the fol- 
lowing integro-differential equation (IDE) for a function / : [0, oofxS* 1 — > K + : 

(22) %(t,e) = -Mf(t,o)+ f f G a {e-e o -v(d i -0 o ))f(t,e i )f(t ) e o )de o de i , 

Ot Jgl J s i 

where M = / gl /(0, 6) dO, a > 0, G a : S 1 -> M+ is the periodic Gaussian with f gi G a (9) d9 = 1, 
and the turning function V : S 1 — > S 1 is odd. 

This IDE describes a jump process in which particles at an old orientation 9 interact over S 1 
with particles in 0i and jump to a new position 8 = 9 + V{Bi — o ). The precision of the jump is 
measured by a. Note that for the IDE V is a turning (therefore it maps to S 1 ), while the function 
V for the TDE is a velocity and takes real values that can be arbitrarily large. If, e.g., V(ip) = ip 
in the IDE, then all solutions converge for t — > oo to the constant solution (Geigant [ ]), while the 
TDE has non-constant stable stationary solutions. 

Both equations preserve mass, positivity, axial symmetry and any periodicity, and both are invari- 
ant under translations and reflections. The S f O(2)-invariance makes linearization and calculation 
of eigenvalues near the stationary homogeneous solution possible, as well as the fast numerical 
calculation of solutions by Fourier transforming the equation into a system of ODEs (see Geigant 
and Stoll [ ] for the IDE). 

Let M = 1. If V = 0, then the solutions of both systems converge to the constant 1 as t — > oo (the 
TDE is the linear diffusion equation, the IDE a linear jump process). If D or a are large compared 
to V, solutions also converge to 1 (Theorem 3.1 for the TDE; Geigant [1] for the IDE). Therefore, 
if V is small, then D and cr, resp., must be very small for pattern formation. On the other hand, 
if V = and D — or a — 0, resp., then |£ = 0, therefore nothing happens. Hence, if V as well 
as D or a are very small, then pattern formation occurs very slowly (if at all). Last but not least, 
if D = or a = but V ^ 0, the limiting equations of both equations have delta distributions as 
solutions. 

This said, we assume that a and V are very small, and we use Taylor expansion in a, V to get 

2 

G a {9 -0 o - V{0i - 9 )) = 6(9 - B ) - V(0i - 9 )5'(9 - 9 ) + °—b"(9 - o ) + 0({a 2 + \V\) 2 ). 
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Plugging this right hand side into (22) yields the transport-diffusion equation (f ) with D = <r 2 /2, 
because 

/ s(e~d )f(e l )f(9 )dd l de = f(9), 

6'(9-9 )(V(9 i -9 )f(9 ))d9 f(9 i )d9 i = f -jL(y{9i- 9 o )f(0 o )) f{9 l )d9 r 

Jgi au e =e 

(V'(di - 8)f(8) - Vid, - 9)f(6))f{9 i ) d9i 

1 

-^{(V*f)(9)f(9)), 



5" (9 - 9 )f{9 )f(9 l ) d0 o d9i = / f(9i) Mi f"(9) = f'{9) 



and 



Different arguments for this derivation are given in Mogilner and Edelstein-Keshet [7] and in Primi 
et al. [9]. 

Similarly, for the eigenvalues Ck of (22) (see Geigant and Stoll [3]), Taylor expansion with small V 
yields 



c fe = -l + 4G CT , fe / cos(7rfcV>) cos(2-Kk{\ij} - V(ip))) dip 
Jo 

w -1 + 4G <T ,fc J cos(7rfcV>) cos(irhl>) dif> + Si^k 2 G a< k / ms| wkr) >\n{ tA-(.-)1 •) <h ' 
= (-1 + G a>k ) + 4n 2 k 2 G^ k ( 2 V{ip) sin(27rfcV) d^ 



°^4ir 2 k 2 / V(ip)sm(2Trkil>)d<ip, 
Jo 

i 

because the fc-th Fourier coefficient G ai k tends to 1 as a — » 0. Because c k — 47rfc J Q 2 V(9) sin(27rfc0) d9 
for D = (see (2)), the signs of the eigenvalues of both models agree for small enough V, D and 
<r (similar arguments were given by 1. Primi, personal communication). We stress again that for 
larger V or D,a, the signs of the eigenvalues may differ. 

But we see for example that for both models there are turning functions V that are negative on 
]0, | [ but lead to non-trivial patterns, see the Example 6.4 in Section 6.3. Especially for the IDE 
this was a surprise to us. Only the eigenvalue of the first mode in the IDE is always negative if V 
is negative, which may correspond to the result of Primi et al. [9] that there are no one-peak like 
solutions for small diffusion if L V{ip) dip is negative somewhere. 

The formulas for the eigenvalues show also that higher modes have larger eigenvalues for the IDE 
[k 2 versus k in TDE). This explains perhaps why in simulations of the IDE at small a we see the 
initial formation of several peak-like maxima much more often than in simulations of the TDE 
with small diffusion D. 

Both equations have limiting equations for D — > and a —> 0, respectively. For a = we have 
G a = Sq, and the limiting equation is 

(23) % (t, 9) = -f(t, 9)+ f f(t, 9 - V(il>))f(t, 9 + ip - V{9)) dyj. 



dt 



,91 
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In Geigant [2] it is shown that for a — > the solutions of the IDE converge to those of the limiting 
equation on fixed finite time intervals. 

For both limiting equations a single peak is a stationary solution, which is linearly stable if V is 
attracting. 'Attraction' in the case of the IDE means < V(9) < 9 for < 9 < ~ (see Geigant [2] 2 ), 
and in the case of the TDE V(9) > 0, V'(0) > for < 9 < 5, see Theorem 5.6. In both equations 
the perturbation may not extend to the opposite side of the peak since particles located there 
cannot turn back (because V{\) =0). 

Two peaks with distance \ whose masses add up to 1 are also a stationary solution for both 
limiting equations because V{\) =0. For the IDE with a — Kang et al. formulate theorems on 
convergence of solutions to two opposite peaks if the initial distribution is sufficiently localized, see 
Theorems 15 and 19 in [(>]. However, there is an estimate in both proofs, namely equations (24) 
and (43), where we cannot follow the argument — they seem to bound a delta distribution by a 
constant. Unfortunately, this estimate is very important for the proofs. We have been informed 
by the authors that an erratum is in preparation. 

For both equations the assumptions for convergence to two opposite peaks are essentially an at- 
tracting shape of V near (V > to the right of 0, V'(0) > 0, and for the IDE additionally 
V'(0) < 1 near 0) and near \ {V < to the left of \, V'{\) > 0, and for the IDE additionally 
V'{\)<1). 

It is important to see that in both limiting equations there is no 'mass selection' toward equal 
masses of the peaks. The open question is then on what time scales equalization of the peaks 
occurs when a or D, resp., are positive. 

The central differences between the two limiting equations for the TDE and IDE are as follows. 

• n > 2 initial peaks, i.e., /(0, •) = Y^k=i m kH- ~ @k) with rrik > 0, do not keep that 
form for the IDE (e.g., starting with two peaks in 9\,9 2l particles jump also to positions 
6>i + V{9 2 - 9f) and 9 2 + V(9 1 - 9 2 )). 

• n > 3 peaks — even if equidistant and with equal masses — are in general not a stationary 
solution for the IDE. 

Therefore, the IDE does not allow the 'peak game' (see Section 5.1; terminology by Mogilner 
et al. [8]). Only if V(^) = for j = 1, . , . ,n — 1, then n equidistant peaks with arbitrary masses 
are a stationary solution of equation (23). It is an educated guess that they are locally stable up 
to redistribution of mass and reorientation if < V 1 (^) < 1 holds for < j < n — 1. 
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